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Abstract 

The models of nn transitions in the medium based on unitary S-matrix are considered. 
The time-dependence and corrections to the models are studied. The lower limits on the 
free-space nn oscillation time are obtained as well. 
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1 Introduction 



In the standard calculations of ab oscillations in the medium [1-3] the interaction of particles a 
and b with the matter is described by the potentials Ua,b (potential model). Imf/f, is responsible 
for loss of 6-particle intensity. In particular, this model is used for the nn transitions in a 
medium [4-10] followed by annihilation: 

n^n^ M, (1) 

here M are the annihilation mesons. 

In [9] it was shown that one-particle model mentioned above does not describe the total ab 
(neutron-antineutron) transition probability as well as the channel corresponding to absorption 
of the 6-particle (antineutron). The effect of final state absorption (annihilation) acts in the 
opposite (wrong) direction, which tends to the additional suppression of the nn transition. The 
yS-matrix should be unitary. 

In [11] we have proposed the model based on the diagram technique which does not contain 
the non-hermitian operators. Subsequently, this calculation was repeated in [12]. However, 
in [13] it was shown that this model is unsuitable: the neutron line entering into the nn 
transition vertex should be the wave function, but not the propagator, as in the model based 
on the diagram technique. For the problem under study this fact is crucial. It leads to the 
cardinal error for the process in nuclei. The nn transitions in the medium and vacuum are not 
reproduced at all. If the neutron binding energy goes to zero, the result diverges (see Eqs. (18) 
and (19) of Ref. [11] or Eqs. (15) and (17) of Ref. [12]). So we abandoned this model [13]. (In 
the recent manuscript [14] the previous calculations [11,12] have been repeated. The model and 
calculation are the same as in [11,12]. Unfortunately, several statements are erroneous [15], in 
particular, the conclusion based on an analogy with the nucleus form-factor at zero momentum 
transfer (for more details, see [15]).) 

In [16] the model which is free of drawbacks given above has been proposed (model b in the 
notations of present paper). However, the consideration was schematic since our concern was 
only with the role of the final state absorption in principle. In Sect. 2 this model as well as the 
model with bare propagator are studied in detail. The corrections to the models (Sect. 3) and 
time-dependence (Sect. 4) are considered as well. In addition, we sum up the present state of 
the investigations of this problem (Sect. 5). 

The basic material is given in Sects. 2 and 5. 



2 



2 Models 



First of all wc consider the antineutron annihilation in the medium. The annihilation amplitude 
Ma is defined as 



</0| Texp(-i J dxn(x)) - 1 |Onp>= A^(27r)^(5^(p/ - Pi)M„. 



(2) 



Here "H is the Hamiltonian of the n-medium interaction, | Oup > is the state of the medium 
containing the n with the 4-momentum p = (e, p); < / | denotes the annihilation products, N 
includes the normalization factors of the wave functions. The antineutron annihilation width 
r is expressed through M^. 

r~ j d^\Ma\^ . (3) 

For the Hamiltonian % we consider the model 

H{t) = j (fxn{x) = Ha{t) + V, (4) 

where Ha is the effective annihilation Hamiltonian in the second quantization representation, 
V is the residual scalar field. The diagrams for the model (4) are shown in Fig. 1. The first 
diagram corresponds to the first order in "H^ and so on. 



n 



n 



n 



M 



n 



n ^ 

-^o-V-o^ - 

' 'v ^ 



n 



/ 

■x 



n 



n 



'V ^ 



Figure 1: Antineutron annihilation in the medium. The annihilation is shown by a circle 



Consider now the process (1). The neutron wave function is 

n{x) = il^^^^ exp{—ipx). 



(5) 



Herep = (e, p) is the neutron 4-momentum; e = p^/2m-|-C/„, where C/„ is the neutron potential. 
The interaction Hamiltonian has the form 



Hi = Hnn + H, 



(6) 



3 



Hnn{t) = j d^x{enn^n{x)'^n{x) + H.C.) (7) 

Here Hnn is the Hamiltonian of nn conversion [6], is a small parameter with = l/r, 
where r is the free-space nn oscillation time; m„ = nin — ra. In the lowest order in H^n the 
amplitude of process (1) is uniquely determined by the Hamiltonian (6): 

M = ennGoMa, (8) 

= e-p2/2m-[/„ + 20' 
Pn = P; = Here Go is the antineutron propagator. The corresponding diagram is 
shown in Fig. 2a. The annihilation amplitude Ma is given by (2). where "H = Tia + V^^n^fj. 
Since Ma contains all the n-medium interactions followed by annihilation including antineutron 
rescattering in the initial state, the antineutron propagator Gq is bare. Once the antineutron 
annihilation amplitude is defined by (2), the expression for the process amplitude (8) rigorously 
follows from (6). For the time being we do not go into the singularity Go ~ 1/0. 

/ / 
- ' n n ' 



Figure 2: a nn transition in the medium followed by annihilation. The antineutron annihilation 
is shown by a circle, b Same as a but the antineutron propagator is dressed (see text) 

One can construct the model with the dressed propagator. We include the scalar field V in 
the antineutron Green function 

G. = G. + GoVG. + ...= |^^^ = -i = 4, (10) 

Tj = V, where S is the antineutron self-energy. The process amplitude is 

M = enuGdMb, (11) 

GdMh = GoMa (see Fig. 2b). The block in the square braces shown in Fig.l corresponds to 
the vertex function M^. The models shown in Figs. 2a and 2b we denote as the models a and 
b, respectively. 
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In both models the interaction Hamiltonians Hj and unperturbed Hamiltonians are the 
same. If S — )■ 0, the model b goes into model a. In this sense the model a is the limiting case 
of the model b. 

We consider the model b. For the process width Tb one obtains 

2 2 

Tb = j d^\M ^iVi jd^\Mb P= ^T', (12) 

T' = Ni J d<l) \ (13) 

where T' is the annihilation width of n calculated through the M5 (and not Ma). The nor- 
malization multiplier A^^i is the same for Fb and F'. The vertex function Mf, is unknown. (We 
recall the antineutron annihilation width F is expressed through the amplitude Ma.) For the 
estimation we put 

Mb = M„, F' = F. (14) 

This is an uncontrollable approximation. 

The time-dependence is determined by the exponential decay law: 

Wb{t) = 1 - e-^*-* ^ Fbt = |f Fi. (15) 

Equation (15) illustrates the result sensitivity to the vahic of parameter E. 

On the other hand, for nn transitions in nuclear matter the standard calculation gives the 
inverse F-dependence [6-9] 

F/2 Ae^-t 
'^••"W-^^«' (ReC/.-C/„W(r/2)^ "-^- 

where Un is the antineutron optical potential. The wrong F-dependence is a direct consequence 
of the inapphcabihty of the model based on optical potential for the calculation of the total 
process probability [9]. (The above-mentioned model describes the probability of finding an 
antineutron only.) 

Comparing with (15), one obtains 

where the values F = 100 MeV and E = KeUn — Un = 10 McV have been used. Strictly 
speaking, the parameter E is uncertain. We have put E = ReC^ — C/„ = 10 MeV only for 
estimation. 

The model b leads to an increase of the nn transition probabihty. The lower hmit on the 
free-space nn oscillation time r ^;^ increases as well: 

T^in = (3.5 - 7.5) • 10« s. (18) 
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This limit exceeds the previous one (see, for example, Refs. [5-8]) by a factor of five. If S ^ 0, 
Wh rises quadratically. So Eq. (18) can be considered as the estimation from below. 

We return to the model shown in Fig. 2a. We use the basis (n, n). The results do not depend 
on the basis. A main part of existing calculations have been done in n — n representation. 
The physics of the problem is in the Hamiltonian. The transition to the basis of stationary 
states is a formal step. It has a sense only in the case of the potential model H — Hp^t — 
ReUfi — Un — ir/2 =const., when the Hamiltonian of n-medium interaction is replaced by the 
effective mass H — > Hp^t — rUes because the Hermitian Hamiltonian of interaction of the 
stationary states with the medium is unknown. Since we work beyond the potential model, the 
procedure of diagonalization of mass matrix is unrelated to our problem. 

The amplitude (8) diverges 

M^ennGoMar^^. (19) 

(See also Eq. (21) of Ref. [13].) These are infrared singularities conditioned by zero momentum 
transfer in the nn transition vertex. (In the model b the effective momentum transfer qq = 
V = takes place.) 

For solving the problem the field-theoretical approach with finite time interval [17] is used. 
It is infrared free. If if = -f^pot) the approach with finite time interval reproduces all the results 
on the particle oscillations, in particular, the nn transition with n in the final state. (Recall 
that our purpose is to describe the process (1) by means of Hermitian Hamiltonian.) 

For the model a the process (1) probability was found to be [10,13] 

Wait) ^ Wf{t) ^ el^f, ri»l, (20) 

where Wf is the free-space nn transition probability. Owing to annihilation channel, Wa is 
practically equal to the free-space nn transition probability. If t — )> oo, Eq. (20) diverges just 
as the modulus (19) squared does. If S — )> 0, Eq. (15) diverges quadratically as well. 

The explanation of the t^-dependence is simple. The process shown in Fig. 2a represents 
two consecutive subprocesses. The speed and probability of the whole process are defined by 
those of the slower subprocess. If l/F <^ t, the annihilation can be considered instantaneous. 
So, the probability of process (1) is defined by the speed of the nn transition: Wa ~ Wf ~ t^. 

Distribution (20) leads to very strong restriction on the free-space nn oscillation time [10,13]: 

= 10^^ yr. (21) 

3 Corrections 

We show that for the nn transition in medium the corrections to the models and additional 
baryon-number-violating processes (see Fig. 3) cannot essentially change the results. First of 
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all we consider the incoherent contribution of the diagrams 3. In Fig. 3a a meson is radiated 
before the nfi transition. The interaction Hamiltonian has the form 



Hj = J + H^- + H. 



(22) 



In the following the background neutron potential is omitted. The neutron wave function is 
given by (5), were p — (po, p) and po — -01, + p^/2m. 
For the process amphtude M30 one obtains 



M^,^gGer,nGM^^~'-\ (23) 

Po - ?o - "T. - (p - q)2/2m + iO' 
where q is the 4-momentum of meson radiated, M^""^) is the amplitude of antineutron annihi- 
lation in the medium in the {n — 1) mesons. As with model a, the antineutron propagator G 
is bare; the n self-energy S = 0. (The same is true for Figs. 3b-3d.) 



n n n 



n n, n n 







n 



n 







n 



m - 

A 











Figure 3: Corrections to the models (a and b) and additional baryon-number- violating processes 
(c and d) 

If g ^ 0, the amphtude M^a increases since G ^ Gg, 

= ^ 2/0 - I (25) 

Po — m — p^/2m 

(The limiting transition g ^ for the diagram 3a is an imaginary procedure because in the 

vertex n — ?> n$ the real meson is escaped and so go > m,^.) The fact that the amplitude 
increases is essential for us because for Fig. 2a g = 0. Due to this Gq ~ 1/0 and Wa ^ Wb. 
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Let and F*^") be the widths corresponding to the Fig. 3a and annihilation width of n in 
the (n) mesons, respectively; F = Y.{n) r*-"^-*- Taking into account that F^"^ is a smooth function 
of y/s and summing over (n), it is easy to get the estimation: 

2 2 

Fsa fti 5 • 10-V%r ^ %F. (26) 

The time-dependence is determined by the exponential decay law: 

Wsa(t) ^ Tsat = %Fi. (27) 

Comparing with (15) we have: Wsa/Wf, — V'^ /m% <^ 1. So for the model b the contribution of 
diagram 3a is negligible. 

For the model a the contribution of diagram 3a is inessential as well. Indeed, using Eqs. 
(27) and (20) we get 

wut) r 



Wait) mir 

where we have put m$ = mTr. Consequently, if 



(28) 



mlt/r » 1, (29) 

and 

rt > 1 (30) 

(see (20)) then the contribution of diagram 3a is negligible. For the nfi transition in nuclei these 
conditions are fulfilled since in this case F ~ 100 MeV and t = Tq = 1.3 yr, where Tq is the 
observation time in proton-decay type experiment [18].) In fact, it is suffice to hold condition 
(30) only because it is more strong. 

In the calculations made above the free-space nn transition operator has been used. This 
is impulse approximation which is employed for nuclear /3 decay, for instance. The simplest 
medium correction to the vertex (or off-diagonal mass, or transition mass) is shown in Fig. 3b. 
In this event the replacement should be made: 

Cnn Cm = e„n(l Ae), (31) 

Ae = esb/^nn, where 635 is the correction to enn produced by the diagram 3b. For the model a 
the limit becomes 

C„ = (l + Ae)10^«yr. (32) 

Obviously, the Ae cannot change the order of magnitude of the n — > n operator 

is essentially zero-range one. The free-space nn transition comes from the exchange of Higgs 
bosons with the mass niH > 10^ GeV [5]. Since niH 3> mw {j^w is the mass of H^-boson), the 
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renormalization effects should not exceed those characteristic of nuclear (3 decay which is less 
than 0.25 [19]. So the medium corrections to the vertex are inessential for us. The same is true 
for the model b. 

Consider now the baryon-number-violating decay n — > n$ [20] shown in Fig. 3c. It leads 
to the same final state, as the processes depicted in Figs. 2, 3a and 3b. Denoting | q |= g, for 
the decay width F^c one obtains 



9o ~ + ^l- '^^'^ parameter e$ corresponding to the vertex n n$ is unknown and so no 
detailed calculation is possible. 

The baryon-number-violating conversion n — >■ A in the medium [20] shown in Fig. 3d 
cannot produce interference, since it contains i^-meson in the final state. For the rest of the 
diagrams the significant interferences are unlikely because the final states in nN annihilation 
are very complicated configurations and persistent phase relations between different amplitudes 
cannot be expected. This qualitative picture is confirmed by our calculations [21] for p-nuclear 
annihilation. It is easy to verify the following statement: if the incoherent contribution of the 
diagrams 3a-3c to the total nuclear annihilation width is taken into account, the lower limit on 
the free-space nn oscillation time Tmin becomes even better. 

To summarise, the contribution of diagrams 3 is inessential for us. 

4 Time-dependence 

The non-trivial circumstance is the quadratic time-dependence in the model a: Wa ~ t'^. The 
heart of the problem is as follows. The processes depicted by the diagrams 2b and 3 are 

described by the exponential decay law. In the first vertex of these diagrams the momentum 
transfer (Figs. 3a-3c), or effective momentum transfer (Figs. 2b, 3d) takes place. The diagram 
2a contains the infrared divergence conditioned by zero momentum transfer in the nn transition 
vertex. This is unremovable perculiarity. This means that the standard S'-matrix approach is 
inapplicable [10,13,17]. In such an event, the other surprises can be expected as well. From this 
standpoint a non-exponential behaviour comes as no surprise to us. It seems natural that for 
non-singular and singular diagrams the functional structure of the results is different, including 
the time-dependence. The opposite situation would be strange. 

The fact that for the processes with q = the S'-matrix problem formulation {oo, — oo) is 
physically incorrect can be seen even from the limiting case H — 0: if Hj — Hnn (see (6)), the 
solution is periodic. It is obtained by means of non-stationary equations of motion and not 




(33) 
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^'-matrix theory. To reproduce the limiting case H ^ 0, i.e. the periodic solution, we have to 
use the approach with finite time interval. 

If the problem is formulated on the interval {t, 0), the decay width F cannot be introduced 
since F = Z^/^i I Sfi{oo, — oo) p /Tq, Tq — )■ oo. This means that the standard calculation 
scheme should be completely revised. (We would hke to emphasise this fact.) The direct 
calculation by means of evolution operator gives the distribution (20). 

Formally, the different time-dependence is due to ^-dependence of amplitudes. We consider 
Eq. (23), for instance. If q decreases, the amplitude M^a increase; in the hmit g — > it is 
singular (see (25)). The point q — corresponds to realistic process shown in Fig. 2a. The 
i^-dependence of this process is the consequence of the zero momentum transfer. 

The more physical explanation of the i^-dependence is as follows. In the Hamiltonian (22) 
corresponding to Fig. 3a we put H — Hnn — 0. Then the virtual decay n — > n$ takes place. 
The first vertex of the diagram 3a dictates the exponential decay law of the overall process 
shown in Fig. 3a. Similarly, in the Hamiltonian (6) corresponding to Fig. 2a, we put H — 0. 
Then the free-space nn transition takes place which is quadratic in time: Wf{t) = e^^^t^. The 
first vertex determines the time-dependence of the whole process at least for small F. We also 
recall that even for proton decay the possibility of non-exponential behaviour is realistic [22-24]. 

5 Discussian and summary 

The sole physical distinction between models a and b is the zero antineutron self-energy in the 
model a; or, similarly, the definition of antineutron annihilation amplitude. However, it leads 
to the fundamentally different results. 

If E — )■ 0, Wb{t) diverges quadratically. This circumstance should be clarified; otherwise 
the model b can be rejected. The calculation in the framework of the model a gives the finite 
result, which justifies our approach from a conceptual point of view and consideration of the 
model a at least as the limiting case. In reality the model a seems quite realistic in itself. 
Indeed, we list the main drawbacks of the model b. 

1) The approximation — Ma is an uncontrollable one. The value of E is uncertain. 
These points are closely related. 

2) The diagram 2b means that the annihilation is turned on upon forming of the self-energy 
part — V (after multiple rescattering of n). This is counter- intuitive since at low energies 
[25,26] 

CTa > 2.5(T„ (34) 

where aa and as are the cross sections of free-space nN annihilation and fiN scattering, re- 
spectively. The inverse picture is in order: in the first stage of the n-medium interaction the 
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annihilation occurs. This is obvious for the nn transitions in the gas. The model a reproduces 
the competition between scattering and annihilation in the intermediate state [27]. 

3) The time-dependence is a more important characteristic of any process. It is common 
knowledge that the ^-dependence of the process probability in the vacuum and medium is 
identical (for example, exponential decay law (15)). In the model a the ^-dependencies in the 
vacuum and medium coincide: Wa ~ t'^ and Wf ~ t'^. The model b gives Wf, ~ t, whereas 
Wf ~ t^. There is no reason known why we have such a fundamental change. 

A) li H — Un, the model a reproduces all the well-known results on particle oscillations [10] 
in contrast to the model b. 

The model a is free of drawbacks given above. The physics of the model is absolutely 
standard. For instance, for the processes shown in Fig. 3 the antineutron propagators are bare 
as well. 

However, there is fundamental problem in the model a: the singularity of the amplitude 
(19). The approach with finite time interval gives the finite result, which justifies the models 
a and b at least in principle. Nevertheless, the time-dependence Wa ~ t"^ and limit (21) seem 
very unusual. The corresponding calculation contains too many new elements. Due to this 

wc view the results of the model a with certain caution. Besides, due to the zero momentum 
transfer in the nn-transition vertex, the model is extremely sensitive to the S. The process 
under study is unstable. The small change of antineutron self-energy E = -> E = V 7^ 0, or, 
similarly, effctive momentum transfer in the nn transition vertex converts the model a to the 
model b: Wa Wh with Wb ^ Wa- (For the processes with non-zero momentum transfer the 
result is little affected by small change of q.) Although we don't see the specific reasons for 
similar scenario, it must not be ruled out. This is a point of great nicety. 

Finally, the values r^j^ = (3.5 — 7.5) ■ 10^ s and r^j^ — 10^^ yr are interpreted as the 
estimations from below (conservative limit) and from above, respectively. Further investigations 
are desirable. 
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